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The drag and diffusion coefficients of heavy quarks propagating through quark gluon plasma 
(QGP) have been evaluated within the ambit of thermal field theory in Hard Thermal Loop (HTL) 
approximations. It is observed that the magnitudes of both the transport coefficients are changed 
significantly from their values obtained when the t channel divergence is shielded simply by Debye 
mass. The implications of this changes in the transport coefficients on the heavy ion phenomenology 
have been discussed. 



PACS numbers: 

I. INTRODUCTION 

The study of the transport coefficients of strongly coupled system is a field of high contemporary interest both 
theoretically and experimentally. In one hand the calculation of the lower bound on the shear viscosity (77) to entropy 
density (s) ratio (77/s) within the frame work of AdS/CFT model [l| has ignited enormous interests among the 
theorists. On the other hand, the experimental study of the r//s for cold atomic systems and QGP and their similarities 
have generated huge interest across various branches of physics (see [|| for a review). In general, the interaction of 
probes with a medium bring out useful information about the nature of the medium. Since the magnitude of the 
i-C ' transport co-efficients are sensitive to the coupling strength hence these quantities can be adapted as useful quantities 
^4 to characterize a medium. In the context of probing QGP expected to be produced in ultra-relativistic heavy ion 
[ collisions at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) energies, we choose the heavy 
D . quarks (HQs), charm and beauty as probes. That is, we would like to extract the drag and diffusion coefficients of the 
• QGP by studying the propagation of HQs through QGP. Selection of HQs as probes has several advantages, such as 
(i) they are produced very early in the collisions and remain extant throughout the evolution of the QGP. As a result, 
the HQs witness the entire evolution of the system. It is expected that the HQ thermalization time is larger than 
the light quarks by a factor m/T where m is the mass of the HQ and T is the temperature. Therefore, the HQs may 
remain out of equilibrium in QGP. (ii) That is the chances of HQs getting thcrmalized in the system is weaker and 
hence do not dictate the bulk properties of QGP. Moreover, the observed transverse momentum suppression (Raa) 
of leptons originating from the decays of heavy flavours produced in nuclear collisions as compared to those produced 
in proton+proton (pp) collisions at the same colliding energy offer us an opportunity to estimate the drag and 
diffusion coefficients of QGP. In the present circumstances we required to describe the motion of the non-equilibrated 
1 HQs in the background of equilibrated system of QGP. An appropriate foundation is provided by the Fokker-Planck 
\ (FP) equation faQ, which reads as follows: 
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where / stands for the momentum-space distribution of the particle (HQs here) undergoing Brownian motion in the 
thermal bath of QGP. The question of HQ thermalization can be addressed by comparing the solution of FP equation 
with the HQs thermal distribution at any given time. Ai and B^ are related to the drag and diffusion co-efficients 
respectively. The interactions of the HQs with the QGP are incorporated in A{ and Bij. That is Ai and Bij can supply 
the information about the nature of the QGP pj-tl9|. The issue of HQ thermalization in QGP can also be addressed 
experimentally by measuring the elliptic flow (1*2) of leptons from the decays of HQs. Therefore, the evaluation of 
the drag and diffusion coefficients of QGP become extremely important. In the present work we attempt to estimate 
both these coefficients by using the techniques of thermal field theory in hard thermal loop (HTL) approximations. 

We will see below that that drag (diffusion) co-efficients are, essentially, momentum (square of the momentum) 
transfer weighted over the squared interaction matrix element (|Af| 2 ). This indicates that an accurate evaluation of 
\M\ 2 is of vital importance. Present work contemplates on determining elastic transport co-efficients of QGP being 
probed by heavy quarks (charm and bottom) within the ambit perturbativc QCD (pQCD) and HTL approximations. 
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The two main elastic process which contribute to the transport coefficients are: Q + q —> Q + q and Q + g —> Q + g. 
Here Q (q) stands for heavy (light) quarks and g denotes gluon. The \M\ 2 for these processes contain ^-channel 
divergence which are normally regulated by introducing thermal mass (to_d) for the exchanged gluons [I2L i.e. 
by replacing t by t — m 2 D in the denominator of the matrix elements. In the present work, instead of shielding the 
divergences simply by (static) Debyc mass we will use the HTL approximated gluon propagator in the t-channel 
diagrams in a self-consistent way. 

The paper is organized as follows. In the next section the general expressions for the drag and diffusion coefficients 
are outlined. In section III we briefly discuss the main steps for evaluating the gluon propagators in HTL approxima- 
tions. Section IV is devoted for presenting results on the drag and diffusion coefficients. The summary and conclusions 
of the present work is presented in section V. The appendix contains the detailed derivation of the matrix elements 
required for the evaluation of drag and diffusion coefficients without small angle approximation. 



II. THE DRAG AND DIFFUSION CO-EFFICIENTS 



In terms of the transition rates the collision integral of the Boltzmann transport equation can written as 









-J 

collisions J 



= J rf 3 k[ w (p + k,k)/(p + k)- w (p,k)/(p)]. 



(2) 



where w(p, k) is collision rate, say for the processes, Q(p) + g(q) — > Q(p — k) + g(q + k). Using Landau approxima- 
tion i.e. by expan ding w ( p + k, k) in powers of k and keeping upto quadratic term, the Boltzmann transport equation 
can be written as fl2, [l7| 



dt dpi 
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where the kernels are defined as 



e? 3 kw(p, k)ki 



(4) 



and 



d 3 kw(p, k)kikj 



(5) 



For I p |— > 0, Aj — > jpi and Bij — > DSij where 7 and D stand for drag and diffusion co-efficients respectively. The 
drag and diffusion coefficients have recently been evaluated within the ambit of AdS/CFT [2(| and pQCD [2l| and 
their importance for jet quenching have been discussed. Eq.|3]is a nonlinear intcgro-diffcrcntial equation known as the 
Landau kinetic equation. The appearance of parton distribution in the expression for uj makes Eq.|3]a non- linear one. 
For the problem under consideration one of the colliding partners (light quarks or gluons) is in equilibrium. In such 
a situation the distribution function which appears in w can be replaced by thermal distribution. As a consequence 
Eqj3] becomes a linear partial differential equation, known as Fokker-Planck (FP) equation. The T dependence of the 
transport coefficients enter through the thermal distribution appearing in lo. 

As mentioned above the drag and diffusion coefficients are related to the expressions for A. t and These 
coefficients can be calculated from the following expression [l2[ with appropriate choice of the function F(p'), 
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where gQ is the HQ degeneracy, F(p' = p — k) is a function of p, q and CM frame scattering angles and cos\ can be 
obtained from, 



s = p 2 + q 2 + 2{E p E q - \p\\q\cosx) 



(7) 



3 



# c .m. an d 4>c.m. arc polar and azimuthal angles of q respectively. Drag (7) can be obtained by the following replacement 
in Eq. M 



For determining diffusion (D) we substitute, 



F(P') = i 



{P-P'f 



(8) 



(9) 



III. HTL APPROXIMATION AND RESUMMED GLUON PROPAGATOR 



As discussed before the calculation of drag and diffusion co-efficients involve the evaluation of amplitudes for 
processes like Q + q — >• Q + q and Q + g — >• Q + g 22] . The amplitudes from bare perturbation theory contains 
t-channel divergence due to low four-momentum. Q = (w, q) gluon exchange. This divergence can be regulated 
by introducing thermal mass of gluon. Here, we study the HTL approximations [23[ and re-summation of gluon 
propagator which will enable us to regulate the t-channel divergence in a self-consistent way and hence will lead to 
comparatively more reliable values of the transport co-efficients. 

Our aim is to find out HTL approximated self-energy of gluon which goes as an input to the resummed gluon 
propagator to be used as effective thermal propagator regularizing the t channel divergence. The gluon self-energy 
in HTL approximation is discussed in detail in Ref. [U, [2f| . In this section we give only an outline of the scheme. 
There are four diagrams which contribute to gluon self-energy (Fig. [1]). The loop integrations can be written down 



K 




(d) 



FIG. 1: Feynman diagrams contributing to gluon self-energy, (a)ghost-gluon loop, (b)four-gluon vertex, (c) quark-antiquark 
pair creation, (d)three-gluon vertex. 

easily if we keep in mind that the loop- momentum, K — (fco, k) is 'hard' compared to external gluon momentum, i.e. 
K » Q which enables us to use simplified three-gluon vertex [24|. Our goal will be to find out T 2 contributions of 
self-energy because the momentum integration is cut-off at the momentum scales T due to the presence of thermal 
distribution function. That is we can take-up the momentum integration J a kdk which blows up as k — > 00, but, 
insertion of distribution function makes it finite even at k — > 00, 

2 T 2 

kf(k)dk = —— (10) 
li 
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The leading contribution in Eq. [10] is given by k ~ T. If we are interested in high-temperature limit we can assume 
|fc| >> \q\ and approximate related quantities accordingly [24 1. 

The effective gluon propagators evaluated in one loop order in HTL approximation enter in the transport coefficients 
evaluated for the processes displayed Figs. [8] and [9] 

For low momentum transfer (i.e.~ gT where g is the colour charge, g = \/47ra7, a s is the strong coupling), one has 
to use the resummed propagator [26j |. The resummed gluon propagator is given by: 

A"" = + *-± + (a - 1) ^ W (11) 

—Q 2 + U T -Q 2 + Hl ' Q 2 1 J 

will need HTL approximated Hl and Ht- The expressions for and in the context of effective propagator are 
discussed in the Appendix. 



IV. RESULTS 



The drag and diffusion coefficients due to elastic collision have been evaluated by using the HTL approximated 
gluon propagators discussed above, they are denoted by jhtl and Dhtl respectively. The corresponding quantities 
evaluated by shielding the divergence in the i-channel gluon propagator by Debye mass will be denoted by 7 and 
D. We use zero-temperature Fermionic propagator for heavy quarks because they remain out of equilibrium in the 
system under consideration. 

The drag co-efficients of charm and beauty quarks moving in QGP has been evaluated by applying HTL resummed 
gluon propagator. The variations of drags with temperature for HQs with momentum, p = 1 GeV as a function 
of temperature are displayed in Fig. [5] The results clearly indicate an enhancement and rapid variation of jhtl 
compared to 7. The increase is more prominent for charm than beauty. We have explicitly checked that in the static 
limit (uj/q = x — > 0) the "/htl approaches 7. The variation of "/htl and 7 with momentum is depicted in Fig. [3] 
for T = 300 MeV. The "/htl is greater than 7 for the entire momentum range considered here. Again, drag being 
the measure of energy loss, increase in drag results in more suppression of heavy flavours measured at RHIC and 
LHC energies. From Fig. [2] we observe that at 400 MeV temperature the jhtl for charm quark is ^33% more 
than the 7. Whereas, the corresponding difference is ~25% for bottom quarks. We also observe that this difference 
increases with the increase in temperature. The momentum dependence of drag is distinctly affected if we consider 
the HTL resummation technique. For a 5 GeV charm the ^htl is ^50% more than 7. For higher momenta (10 
GeV) the difference reduces to ^45%. In Fig. @] the interaction rate of the HQs with the QGP is shown as function 
of temperature. The rate is increased when the HTL corrections to the t-channel gluon is taken into account. The 
inclusion of effective thermal gluon propagator increases the likelihood of charm or bottom being equilibrated with 
the medium. Therefore, this result will have important implications for elliptic flow of heavy flavours. 

In Figs. [5] and [6]the diffusion coefficients Dhtl and D are plotted with T (for p = 1 GeV) and p (for T — 300 MeV) 
respectively. Similar to drag we observe that the difference between Dhtl and D at high momentum is significant. 
Diffusion co-efficients seem to be more sensitive to the use of effective propagator in a sense that we observe ^100% 
change between the Dhtl and D at T = 400 MeV and this difference increases with T. Though unlike drag, this 
difference is not much (~3.5%) for a difference in charm and beauty quark masses. The momentum dependence of 
diffusion is significantly modified too. A 5 GeV charm diffuses ^80% more when the exchanged gluon contains the 
HTL approximated changes in its spectral function. These changes in drag and diffusion co-efficients originate from 
the spectral modification of the t-channel gluons due to its interaction with the thermal bath. In the static limit 
7rr — > and 7T£ — > m 2 D . The appearance of non-zero ttt makes jhtl larger than 7. 

Now some comments on the magnitude of the values of Dhtl and D are in order here. The value of the diffusion 
coefficients in spatial co-ordinate, D^ TL can be estimated from the value of drag by using the relation D^ TL = 
T {rwy htl)~^ ■ In Fig.[7]wc plot D^ TL multiplied by the inverse of the thermal dc Broglic length, A = l/(2nT). The 
results clearly indicate that the D x remains well above the quantum bound. 



V. SUMMARY AND CONCLUSION: 



In summary, we have taken into account the thermal modifications of the gluon spectral function obtained within 
the ambit of HTL approximations in evaluating the drag and diffusion coefficients of HQs propagating through the 
QGP. The deviations between jhtl and 7 and Dhtl and D is found to be substantial. The enhanced drag will 
result in higher suppression of the HQ momentum spectrum. The increase in drag will also enhance the chances of 
HQ getting equilibrated with the bulk of the system. These results will have crucial consequences on the observables 
like nuclear suppressions, Raa{pt) and elliptic flow, vi of heavy flavours measured at RHIC and LHC experiments. 
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FIG. 2: (Color online) Variation of drag of heavy quarks of momentum 1 GeV with temperature. 
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FIG. 3: (Color online) Variation of drag of heavy quarks with momentum in a QGP bath of temperature 300 MeV. 

VI. APPENDIX: 

In this appendix we demonstrate the calculation of matrix elements of the processes Qq — > Qq (Fig. by applying 
HTL approximation. The matrix elements for the process, Qg — > Qg (Fig. [9j> can be evaluated in a similar way. 
Before decomposing the gluon self energy into longitudinal and transverse components in the medium we define the 
following useful quantities J24[. Let m m be the fluid four- velocity (with normalization condition u^u^ = 1) then any 
four- vector Q**, can be decomposed into a part which is parallel to fluid velocity and the other perpendicular to fluid 
velocity as: 



Q, 



U! = Q.u 
- u^(Q.u) 



(12) 



where 



Q 2 = u? - q 2 
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FIG. 4: (Color online) Temperature variation of the interaction rate of heavy quarks with QGP. 
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FIG. 5: (Color online) Variation of diffusion of heavy quarks of momentum 1 GeV with temperature. 



Q 2 = -q 2 



(13) 



Eqs. H121 and [TBI are valid in the local rest frame of fluid, i.e. in a frame where u = (1, 0). Similarly a tensor orthogonal 
to can be defined, 



(14) 



The full gluon propagator with momentum Q is obtained from the vacuum polarization by using Dyson- Schwinger 
equation: 



-Q 2 + U T -Q 2 + U L 

T 



+ {a-iy w 



Q 2 



(15) 



where a is a gauge-fixing parameter. The longitudinal tensor V^ v and the transverse tensor VT' are defined as [2f 



■p^ v = - 

L D2„2 



(wQ'*-QV)(wQ , '-QV) 



(16) 
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FIG. 6: (Color online) Variation of diffusion of heavy quarks with momentum in a QGP bath of temperature 300 MeV. 
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FIG. 7: Comparison of co-ordinate space diffusion of charm quark with thermal de Broglie wavelength. 



which arc orthogonal to as well as to each other, i.e. 



But, 



V? p V ivp = Vt ,i = L/T 
The free gluon propagator at zero temperature is 

Q»Q V \ 1 



V Q 2 J Q 2 

The transverse and longitudinal self-energies are given by 

n L (Q) = (l-x 2 )w L (x), Tl T (Q) = ir T (x) 



(17) 

(18) 
(19) 

(20) 

(21) 
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where x = oo/q and scaled self-energies ttt and ttl are given by [2j 



tt t (x) = m 2 D 



X X , 2 x > I *^ 



" 2 \ 

i— 2(1 — x ) 
4 y 



n L (x) = m 



D 



x 1 + a; ,7T 



(22) 



Non-zero imaginary parts of the self-energies corresponds to different processes taking place in the medium. 




FIG. 8: Qq — > Qq Feynman diagram. Bold lines are for heavy quarks(Q). 
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FIG. 9: Qg — > Feynman diagrams 



A. Qg — > Qq Matrix Element from HTL approximation: 

The ^-channel gluons that appears in the matrix element for the process Qq — > Qq (Fig[8]) is taken as the effective 
gluon propagator obtained by HTL approximation. In Figs. [5] and [5] the effective propagator will be denoted by a 
solid circle. We can write the amplitude in Feynman Gauge(a = 1) from fig. |8]as, 

- iM t = «(p3)(-t57 /1 i&)u(pi) H A ^] 

u(P4)(^g^t? k )u( P2 ) (23) 

where g is strong coupling and g 2 = 4ira s . k,l (i ^ j, k ^ I) are quark colours and 'a' is the colour of intermediary 
gluon with polarizations fj,, v. After squaring and averaging over spin and colour and using Eq. [15] we get, 
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\Mq£ 
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(24) 



where Cq 9 = | is the color factor, Q 2 = t = (pi 
following relations. 



P3j 



A = < 2 — t(i?e.n T + i?ell i ) + i?e.n T n^ and we have used the 



A""A*; = i ^ + ^ — ? 

p (t-u T f (t-n L f 



IAI 



a"" a: 



(t-n T ) 2 (t-U L f 



(25) 



Using cqs. [16j H2]we can show that 



P\-Vl-P2 = P3-Vl-P4 = P4-Vl-P\ = P2-V L -P3 (26) 
where all the calculations are done in the rest frame of fluid element. 
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